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Abstract 

We have developed a new methodology and a time-domain software package for the estimation 
of the oscillation frequency and the phase noise spectrum of non-linear noisy microwave circuits 
based on the direct integration of the system of stochastic differential equations representing the 
circuit. Our theoretical evaluations can be used in order to make detailed comparisons with the 
experimental measurements of phase noise spectra in selected oscillating circuits. 
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I. INTRODUCTION 



Electronically tuned microwave oscillators are key components used in a wide variety of 
microwave communications systems [lj]. The phase of the output signal exhibits fluctuations 
in time about the steady state oscillations giving rise to Phase Noise a very important char- 
acteristic that influences the overall performance especially at higher microwave frequencies. 

In order to understand the oscillator phase behaviour, a statistical model for a non-linear 
oscillating circuit has to be developed and presently, no accurate theoretical model for phase 
noise characterization is available because of the particularly difficult nature of this problem. 
This is due to the hybrid nature of non-linear microwave oscillator circuits where distributed 
elements (pertaining usually to the associated feeding or resonator circuits) and non-linear 
elements (pertaining usually to the amplifiying circuit) have to be dealt with simultaneously 

Q. 

The main aim of this report is to establish a theoretical framework for dealing with the 
noise sources and non- linearities present in these oscillators, introduce a new methodol- 
ogy to calculate the resonance frequency and evaluate the time responses (waveforms) for 
various voltages and currents in the circuit without or with the noise present. Once this 
is established, the phase noise spectrum is determined and afterwards the validity range 
of the model is experimentally gauged with the use of different types of microwave oscilla- 
tors . This report is organised in the following way: Section II covers the theoretical 
analysis for the oscillating circuit, reviews noise source models and earlier approches. Sec- 
tion III presents results of the theoretical analysis and highlights the determination of the 
resonance frequency for some oscillator circuits without noise. In section IV, phase noise 
spectra are determined for several oscillator circuits and section V contains the experimen- 
tal results. The Appendix contains circuit diagrams and corresponding state equations for 
several non-linear oscillator circuits. 

II. THEORETICAL ANALYSIS 

In standard microwave analysis, it is difficult to deal with distributed elements in the time 
domain and difficult to deal with non-linear elements in the frequency domain. Non- linear 
microwave oscillator circuits have simultaneously non- linear elements in the amplifying part 
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and distributed elements in the resonating part [Non-linearity is needed since it is well known 
that only non-linear circuits have stable oscillations]. 

Before we tackle, in detail, the determination of the phase noise, let us describe the 
standard procedure for dealing with the determination of resonance frequency of non-linear 
oscillator circuits: 

• The first step is to develop a circuit model for the oscillator device and the tuning 
elements. The equivalent circuit should contain inherently noiseless elements and noise 
sources that can be added at will in various parts of the circuit. This separation is 
useful for pinpointing later on the precise noise source location and its origin |2]. The 
resulting circuit is described by a set of coupled non- linear differential equations that 
have to be written in a way such that a linear sub-circuit (usually the resonating part) 
is coupled to another non-linear sub-circuit (usually the oscillating part). 

• The determination of the periodic response of the non- linear circuit. 

• The third step entails performing small signal ac analysis (linearization procedure) 
around the operating point. The result of the ac analysis is a system matrix which is 
ill-conditioned since a large discrepency of frequencies are present simultaneously (One 
has a factor of one million in going from kHz to GHz frequencies) . The eigenvalues of 
this matrix have to be calculated with extra care due to the sensitivity of the matrix 
elements to any numerical roundoff jjj. 

We differ from the above analysis, by integrating the state equations directly with 
standard/non- standard Runge-Kutta methods adapted to the non- stiff /stiff system of or- 
dinary differential equations. The resonance frequency is evaluated directly from the wave- 
forms and the noise is included at various points in the circuit as Johnson or Shot noise. 

This allows us to deal exclusively with time domain methods for the noiseless/noisy non- 
linear elements as well as the distributed elements. The latter are dealt with through an 
equivalence to lumped elements at a particular frequency. 

As far as point 3 is concerned, the linearization procedure method is valid only for small- 
signal analysis whereas in this situation, we are dealing with the large signal case. 

Previously, several methods have been developed in order to find the periodic response. 
The most well established methods are the Harmonic balance and the piecewise Harmonic 



balance methods j(3]. Schwab ^| has combined the time-domain (for the non-linear amplifier 
part) with the frequency domain (for the linear resonating part) methods and transformed 
the system of equations into a boundary value problem that yields the periodic response of 
the system. 

III. TIME RESPONSES OF NON-LINEAR OSCILLATORS AND RESONANCE 
FREQUENCY DETERMINATION 

For illustration and validation of the method we solve 6 different oscillator circuits (The 
Appendix contains the circuit diagrams and the corresponding state equations): 

• The standard Van der Pol oscillator. 

• The amplitude controlled Van der Pol oscillator. 

• The Clapp oscillator. 

• The Colpitts oscillator. 

• Model I oscillator. 

• Model II oscillator. 

We display the time responses (waveforms) for various voltages and currents in the at- 
tached figures for each of the six oscillators. All oscillators reach periodic steady state almost 
instantly except the amplitude controlled Van der Pol (ACVDP) and the Colpitts circuits. 
For instance, we need, typically, several thousand time steps to drive the ACVDP circuit 
into the oscillatory steady state whereas several hundred thousand steps are required for the 
Colpitts circuit. Typically, the rest of the circuits studied reached the periodic steady state 
in only less a couple of hundred steps. 

IV. PHASE NOISE SPECTRUM EVALUATION 

Once the oscillating frequency is obtained, device noise is turned on and its effect on the 
oscillator phase noise is evaluated. All the above analysis is performed with time domain 
simulation techniques. 



Finally, Fourier analysis is applied to the waveform obtained in order to extract the power 
spectrum as a function of frequency. Very long simulation times (on the order of several 
hundred thousand cycles) are needed since one expects inverse power-law dependencies on 
the frequency 

We use a special Stochastic time integration method namely the 2S-20-2G Runge-Kutta 
method developed by Klauder and Peterson, and we calculate the PSD (Power Spectral 
Density) from the time series obtained. 

It is worth mentioning that our methodology is valid for any type of oscillator circuit and 
for any type of noise (Additive White as it is in Johnson noise of resistors, Mutiplicative 
and Colored or l/f a with a arbitrary as it is for Shot noise stemming from junctions or 
imperfections inside the device). In addition, the approach we develop is independent of the 
magnitude of the noise. Regardless of the noise intensity we evaluate the time response and 
later on the power spectrum without performing any perturbative development whatsoever. 
Recently, Kartner (jf developed a perturbative approach to evaluate the power spectrum 
without having to integrate the state equations. His approach is valid for weak noise only 
and is based on an analytical expression for the power spectrum. Nevertheless one needs to 
evaluate numerically one Fourier coefficient g^o the spectrum depends on. 

V. EXPERIMENTAL VERIFICATION 

Microwave oscillators are realised using a very wide variety of circuit configurations 

and resonators. We plan to design, fabricate and test microstrip oscillators with GaAs 

n 

MESFET devices with coupled lines and ring resonators (5j. The measured phase noise of 
these oscillators will be compared with the theoretical prediction from the above analysis. 
We also plan to apply the above analysis to the experimental phase results obtained from 



various e 
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ectronically tuned oscillators that have been already published in the literature 
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two additional circuits (Model I and II) to test the software. 
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APPENDIX 

A. Van Der Pol oscillator 

State-Space Equations of Van der Pol oscillator: 
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di L Vc 
dt ~ L 

= g(-nV c i L 2 -i L + fiV c ) 



Define: 



x 1 


Def. 




x 2 







Rewrite Equations 1 and 2 in state-space form: 

Xi = X 2 

X 2 = fiX 2 (l - Xf) - x 1 
B. Amplitude controlled Van Der Pol oscillator 

State-Space equations of Amplitude controlled Van der Pol oscillator 



= - — (i L + Gu-G a (z)u-i r ) 

diL u 
~dt ~ ~L 

dz 1 r . u N 2 , 

Tt = T [i vJ ~ z] 

where G a (z) — G(l — c(z — 1)) is the voltage controlled conductance. 



Choosing the variables: 







U 


x 2 


Def. 


G 






Z 



with the definitions Qo = and i rn = 7 = 7^ where uj is the oscillation 

frequency, we have: 



x 3 



cX 1 (X 3 -l)+X 2 

Q0X1 

- 7 (X 3 - Xl) 



+ It 



(10) 

(11) 

(12) 



C. Clapp oscillator 

The parameters used are: 



r s = on 


R E = 500Q 


r a = moon 


R L = lOOOfi 


C A = lOOpF 


C E = WOpF 


C TE = 3.5nF 


C TA = 7nF 



In addition, we have: C s = 25 pF and L s — 1 mH. 
State-Space Equations of Clapp oscillator: 



dV C E 




dt 




dVcTE 


1 


dt 


Cte 


dVcTA 


1 


dt 


Cta 


dV CA 




dt 




dip 


= jp 


~dt 



■ VcTE — Vqe 



)] 



Re 

Vote — Vce 
Re 

Ra 



)] 



• VcTA — V C A ^ _ V C A 1 



Ra 



R, 



(13) 
(14) 
(15) 
(16) 

(17) 



whereas: 



dt 



R s j P ,Vce 

T V 



L s ReCte 



VcTE^ ,VcTA 



CA^ 



L s RaCta 



where: 



and: 



. 1 + _J_ + _J_ 

L s Cta L s C s L s Cxe 



CE 



i e = A 1 (e A * -1) 



-B 1 B 2 
—BiVca 
B\B 2 



- V C a < ~B 2 or V C a > B 2 

— B 2 < Vca < B 2 
- Vca > B 2 or V CA < -B 2 



D. Colpitts oscillator 



Ri = 350Q 


R 2 = llOJfeft 


R L = 500Q 


L x = lOmH 


L 2 = 30nH 


d = lOpF 


C 2 = 940pF 


C 3 = 2.7nF 


C A = l.hnF 



For the transistor: 



I s = 2.7 10- 16 A 


fa = 5.5 


Pn = 140 


U A = 15V 


U B = 4.3 V 


V = 12V 



State-Space Equations of Coplitts oscillator: 



~dl 
di 2 
~dl 
dV 1 

~dt 

(TV 
dt 

dV3 
dt 

dV A 

~dt 



= ^(Vo-Vi-Rih) 

= Uv,-v 2 ) 

= -^rbi - 12 - *c + {—^ — )] 

l\. • , , Vo-V2-V 3 y 
= -rrfa-iB + i 5 )] 

t/2 ti 2 



1 r . ,V -V 2 -V 3 



[-IB + (- 



Jj. ( 



c 3 

V-V 4 



Ri 



)] 



RlC a 
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Where the transistor currents are given by: 



% c — % ce — t>bc 



IB = iBE + iBC 



I s U BE U BC 

icE = 7nrr(e " T -e* T ) 



(Qb_\ 

I- , U BE 



iBE = -5-(e u t - 1) 

PN 

I s , U BC 

iBC = ~r(e - 1) 
Pi 

Moreover, we have the additional relations: 



ube = V 2 + V 3 

u BC = V 2 + V 3 -V 1 

Qb _ 1 ube_ ubc 
Qw u b u a 



V2+V3 V2+V3-V1 
(e u t — e u t ) 

ic = s d I V2+V3 I EEEIEEK T 

V U b ' U a > 

- Trie »t - 1) 



i B = — (e »t - 1) + — (e »t 
Pat Pj 



Define: 







i\ 


X 2 




i2 


X 3 


Def. 


v 1 


x 4 




v 2 


x 5 




V 3 


x Q _ 




v 4 _ 
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Rewrite Equations in state-space form: 



X x = ^(V - X 3 - R 1 X 1 ) 
X 2 = -j—(X 3 — X4) 

J^2 



I I T x 4 +x 5 -x- i 

X 3 = -[X 1 -X 2 + — (X 6 - X 3 ) + -L (e—^- - 1) 
Oi Kl pi 



I s x 4,+ x b X 4 +X 5 -X 3 



a I X4+X5 _|_ X4+X5—X3 S 
Ub Ua 




I I I x 4 +x 5 

X 5 = — [— (Vb - X 4 - X 5 ) - ^-(e— - 1) 

^3 -n-2 Pat 
J s *4+* 5 -*a 
- 7r(e ur - 1)] 
Pi 



x e - D n {X 3 - X 6 ) 
rtLL/4 



E. Model I 



V = W 


R 1 = 220hQ 


R 2 = lOOOfi 


R 3 = 220ktt 


i? 4 = 2fi 


L = lOfiH 


d = 0.47/iF 


C 2 = 200pF 


C 3 = 200pF 



For the transistor: 
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I a = 2.7 10~ 16 A 


U A = 15V 


U B = 4.3 V 


Gain 


Pi = 5.5 


p N = 140 



State-Space Equations of Model I oscillator: 



di4 
~dt 

dt 



-(V 6 + V 7 -R 4 i 4 -V 5 ) 

ii 

Ci 



dVe = 1 [( V -V 6 -V 7 ) _ { Vs + V7 ) 
dt C*2 Ri R3 



m 1 V -V 6 -V 7 
dt C 3 [[ i?i ' 

- « 4 - IB + IE ~ -5- J 

-^2 



V + V 



The currents and transistor voltages «be and ubc are given by: 

^ = ib + ic 
ube = Vq 
ubc = V 6 + V 7 -V 

Therefore, %b and ic are given by: 



Is , ^ 

Pn 



l) + f(e -r -i; 
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^ ^ U b ^ Ua > 

I . V6+V7-VD 

- 7r(e UT - 1) 
Pi 



Define: 



x 1 




u 


x 2 


Def. 


v 3 


x 3 










v 7 _ 



Rewrite equations in state-space form: 



x 2 



L 

Xi 

Ci 



(X 3 + X 4 — X 2 — R4X1) 




X 4 = ^-[^-(V -X 3 - X 4 ) - X 1 - -^(X 3 + X 4 ) 

03 ill K 3 

X 4 I s 

I s £3 x 3 +x 4 -v 

+ (1 I x 3 , x 3 +x 4 -VbT ( e " T ~ e " T )1 



F. Model II 



i?i = lOOOfi 


R 2 = 82k£l 


R 3 = 680Q 


L = lOOnH 


Ci = 33pF 


C 2 = 33pF 


C 3 = lOpF 


C A = 0.1/iF 
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For the transistor: 



I a = 2.7 10~ 16 A 


U A = 15V 


U B = 4.3 V 


Gain 


p! = 5.5 


(3 N = 140 



State-Space Equations of Model II oscillator: 

The transistor voltages and are given by: 



dV 
~dt 



1 ,, Vo-Vi - 
<V l R 1 



V x 



i2 3 



Vi-V 2 



B.2 

)-ie] 



dV2 1 Vi- 
dt C 2 [[ R 2 



Yi ) + ( Yl 



Vo-V, 



+ ie~ iB\ 



C 3 

±(V 1 -V 2 -V b ) 
1 ( Vi-V2-V 5 ) 

The currents and transistor voltages Ube and Ubc are given by: 

ube = V 2 
ubc = V 2 -Vi 



dt 
die 

dt 
dVs 

dt 



Therefore, %b and %c are given by: 



I S , %L X I S / V2 ~ Vl 

i B = — (e»T - 1) + — (e «r - 1) 

Pn Pi 
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^ ^ u b T 



I s , V 2- V l 

_ ( e u T 
Pi 



-1) 



Define: 







Vl" 


x 2 




^2 


x 3 


Def. 


^3 


x 4 






. X 5. 







Rewrite equations in state-space form: 





X 3 
X 5 



C 3 



-\X\ — X2 — X3) 



1 



C4-R3 



.( Xl -X 2 -X 5 ) 



FIGURE CAPTIONS 
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FIG. 1: Van der Pol oscillator circuit 




-G a (z) 




FIG. 2: Amplitude controlled Van der Pol oscillator circuit 




FIG. 3: Clapp oscillator circuit 
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FIG. 4: Colpitts oscillator circuit 




FIG. 5: Model 1 oscillator circuit 
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I v o = 15 V 



V '1 



v 6 



V .4 



L 



C 3 



'6 
► 



5V < 



v 5 



V 7 Vic 



1 ' >2 +i 3 



V '7 



FIG. 6: Model II oscillator circuit 
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